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►  The  uniform  charging  cell  voltage  curves  (CCVC)  hypothesis  is  proposed. 

►  Cell  capacities  can  be  estimated  by  overlapping  CCVCs  using  CCVC  transformation. 

►  A  simplified  approach  using  voltage-capacity  rate  curve  (VCRC)  is  proposed. 
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►  The  battery  pack  capacity  can  be  precisely  estimated  with  the  proposed  method. 
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Because  of  the  diversiform  driving  conditions  and  the  cell  variations,  it  is  difficult  to  accurately  determine 
battery  pack  capacities  in  electric  vehicles  (EVs)  by  model  prediction  or  direct  measurement.  This  paper 
studies  the  charging  cell  voltage  curves  (CCVC)  for  the  estimation  of  the  LiFeP04  battery  pack  capacities 
in  EVs.  We  propose  the  uniform  CCVC  hypothesis  and  estimate  cell  capacities  by  overlapping  CCVCs  using 
CCVC  transformation.  CCVCs  of  two  LiFeP04  cells  with  large  capacity  difference  are  used  to  verify  the 
hypothesis.  We  further  develop  an  equivalent  simplified  approach  using  voltage-capacity  rate  curve 
(VCRC)  and  implement  Genetic  Algorithm  (GA)  to  find  the  optimum  transformation  parameter  for 
overlapping  VCRCs.  A  small  battery  pack  with  four  LiFeP04  cells  in  series  is  employed  to  verify  the 
method  and  the  result  shows  that  the  estimation  errors  of  both  pack  capacity  and  cell  capacities  are  less 
than  1%.  With  the  proposed  method,  the  battery  pack  capacity  can  be  precisely  estimated  which  could  be 
used  for  the  driving  range  prediction.  Meanwhile,  the  estimated  cell  capacities  in  battery  packs  will 
significantly  support  the  study  of  cell  degradation  and  cell  variations  in  vehicle  driving  conditions. 

Crown  Copyright  ©  2012  Published  by  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

With  rapid  progress  in  power  density  and  energy  density 
increase  of  the  lithium-ion  batteries  over  the  past  decade,  electric 
vehicles  (EVs)  with  lithium-ion  battery  packs  as  power  sources 
have  been  intensively  concerned.  However,  as  one  single  cell 
provides  insufficient  voltage  and  capacity,  it  is  inevitable  to 
construct  battery  packs  with  hundreds  of  single  cells  connected  in 
parallel  and  series  to  meet  the  power  and  energy  required  in  EVs 
[1,2].  If  there  were  no  differences  between  the  cells,  we  could 
simply  consider  the  battery  pack  as  one  cell  with  high  voltage  and 
large  capacity.  Current  researches  on  the  vehicle  level  usually  adopt 
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this  simplified  method  [3-5],  Unfortunately,  due  to  the  inconsis¬ 
tent  manufacturing  process  and  the  inhomogeneous  operating 
environment,  cells  always  have  variations  which  cannot  be  elimi¬ 
nated  [6,7], 

Capacity  variations  of  series-connected  cells  directly  influence 
the  pack  capacity.  The  pack  energy  density  and  the  driving  range  of 
an  EV  depend  on  the  pack  capacity  principally.  Thus,  it  is  of  great 
significant  to  estimate  the  precise  pack  capacity  considering  cell 
capacity  variations.  The  initial  pack  capacity  can  be  actually  calcu¬ 
lated  through  the  experiment  which  is  a  standard  capacity  test 
method  recommended  by  the  manufactory.  Nevertheless,  pack 
capacity  changes  as  cell  capacities  fade.  As  a  result,  the  measured 
initial  pack  capacity  is  meaningless.  For  electric  vehicles,  depleting 
the  remaining  pack  capacity  would  be  good  for  measuring  the  pack 
capacity  during  operating  but  it  would  be  a  disaster  to  the  driver  as 
nobody  is  willing  to  pull  his  car  to  the  charging  station.  Research 
and  literature  into  single  cell  capacity  fade  prediction  are  still  the 
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most  concerned  topics  [8—11],  E.V.  Thomas  et  al.  [8]  and  M.  Ecker 
et  al.  [9]  used  experimental  results  of  several  cells  in  different 
conditions  to  develop  the  lifetime  model.  In  their  experiments,  fade 
rates  of  cell  capacities  were  different  even  when  cells  were  in  the 
same  test  conditions.  Therefore,  it  would  be  very  difficult  and 
inaccurate  to  predict  the  pack  capacity  fade  using  similar  methods 
that  predict  single  cell  capacity.  As  a  consequence,  little  literature 
explores  prediction  methods  for  pack  capacities. 

Because  it  is  almost  impossible  to  measure  the  pack  capacity  in  EVs 
due  to  the  driving  condition  and  also  predicting  the  pack  capacity  is 
inaccurate,  precise  pack  capacity  estimation  in  EVs  is  a  challenging 
task.  Several  cell  capacity  estimation  approaches  were  studied  in 
literature  [12-15],  G.  L.  Plett  [12]  and  H.F.  Dai  et  al.  [13]  used  dual 
Kalman  filtering  for  state  of  charge  (SOC)  and  capacity  estimation.  C. 
Hu  et  al.  [14]  proposed  an  improved  multiscale  framework  with  Kal¬ 
man  filtering  for  SOC  and  capacity  estimation.  The  estimation  results 
were  satisfactory  as  accurate  battery  models  for  cells  were  imple¬ 
mented.  However,  it  is  difficult  to  construct  an  accurate  model  con¬ 
cerning  inconsistent  cells  in  the  battery  pack.  G.  L  Plett  [15]  utilized 
a  total  least  squares  estimation  method  to  estimate  cell  capacity,  and 
no  accurate  battery  models  were  involved  using  this  method.  K. 
Honkura  et  al.  [10]  analyzed  discharge  curves  to  predict  capacity  fade 
of  lithium-ion  batteries.  Nevertheless,  accurate  pack  capacity  esti¬ 
mation  regarding  cell  variations  was  not  focused  in  above  studies. 

This  paper  first  introduces  the  definition  of  battery  “pack 
capacity”  concerning  cell  variations  in  Section  2.  And  subsequently 
we  associate  pack  capacity  with  cell  capacities.  We  propose  the 
uniform  charging  cell  voltage  curve  (CCVC)  hypothesis  in  Section  3: 
when  internal  resistances,  initial  remaining  cell  capacities  and  total 
cell  capacities  of  cells  in  one  batch  are  the  same,  CCVCs  of  the  cells  in 
the  pack  are  the  same,  which  means  that  voltage— capacity  curves 
are  completely  overlapped.  And  we  deduce  that  by  CCVC  trans¬ 
formation,  CCVCs  of  cells  with  different  internal  resistances,  initial 
remaining  cell  capacities  and  total  cell  capacities  can  be  overlapped 
and  cell  capacities  can  be  calculated.  CCVCs  of  two  LiFeP04  cells  with 
large  capacity  difference  are  used  to  verify  the  hypothesis  and  the 
result  shows  that  the  calculated  cell  capacity  is  satisfactory.  We 
further  develop  an  equivalent  simplified  approach  using  voltage- 
capacity  rate  curve  (VCRC)  transformation.  In  Section  4,  Genetic 
Algorithm  (GA)  is  implemented  to  find  the  optimum  transformation 
parameter.  A  small  battery  pack  with  four  LiFePC>4  cells  in  series  is 
employed  to  verify  the  method  and  the  result  shows  that  the  esti¬ 
mation  errors  of  both  pack  capacity  and  cell  capacities  are  qualified. 
With  the  proposed  method,  data  of  CCVCs  can  be  used  to  estimate 
pack  capacities  in  EVs,  which  is  benefit  to  accurate  driving  range 
estimation.  In  addition,  cell  capacities  calculated  in  battery  packs  will 
support  the  study  of  cell  degradation  and  variations  in  vehicle 
driving  conditions,  and  the  calculated  cell  capacities  are  also  the  key 
information  for  cell  equalization  strategies. 


2.  Pack  capacity  and  charging  cell  voltage  curves 

2.1.  Pack  capacity 

Pack  capacity  is  total  available  ampere-hours  released  from 
a  fully  charged  state  of  one  cell  in  the  battery  pack  to  a  fully 


discharged  state  of  one  or  other  cell  at  a  constant  discharge  rate 
of  1/3C  at  room  temperature  of  25  °C.  Pack  capacity  is  commonly 
defined  by  discharging  capacity,  but  regarding  high  coulombic 
efficiency  of  commercial  LiFePC>4  cells  in  EVs,  difference  of 
charging  capacity  and  discharging  capacity  is  negligible.  As  dis¬ 
charging  process  depends  on  EV’s  driving  condition  which  is 
volatile,  while  the  charging  process  is  fixed  at  a  designed 
charging  current,  pack  charging  capacity  as  pack  capacity  is 
satisfactory  in  EVs. 

When  one  or  some  cell  is  fully  discharged  (marked  as  Cell  A) 
during  discharging  without  cell  equalization,  other  cells  still  have 
some  remaining  capacities  because  of  the  difference  in  cell  capacity 
and  SOC.  But  the  battery  pack  is  not  allowed  to  discharge;  other¬ 
wise  it  would  greatly  reduce  the  life  of  the  battery  pack,  or  even 
raise  severe  safety  problems.  Similarly,  when  one  or  some  cell  is 
fully  charged  (marked  as  Cell  B,  and  Cell  B  may  be  Cell  A)  during 
charging,  other  cells  are  still  not  fully  charged,  but  the  battery  pack 
is  not  allowed  to  charge.  The  charged  capacity  during  charging 
process  from  fully  discharged  Cell  A  to  fully  charged  Cell  B  can  be 
called  pack  capacity  as  defined.  When  Cell  B  is  Cell  A,  pack  capacity 
is  exactly  the  capacity  of  Cell  A.  The  releasable  capacity  of  a  battery 
pack  is  the  minimum  remaining  cell  capacity  and  the  absorbable 
capacity  is  the  minimum  capacity  that  can  be  charged.  Pack 
capacity  then  can  be  considered  as  the  sum  of  the  minimum 
remaining  cell  capacity  and  the  minimum  capacity  that  can  be 
charged.  This  can  be  formulated  as 


CpackM  =  min(Cr(t))  +  min((l  -  SOC(t))  -C(t))  (1) 

where  Cpack(t)  is  the  pack  capacity  at  its  lifetime  t,  Cr(t)  is  the 
remaining  cell  capacity,  SOC(t)  is  the  cell  current  state  of  charge  and 
C(t)  is  the  cell  capacity.  When  the  cell  with  the  minimum  remaining 
cell  capacity  has  the  minimum  capacity  that  can  be  charged,  the 
pack  capacity  is  equal  to  the  cell  capacity,  and  the  capacity  of  the 
cell  is  also  the  minimum  in  the  battery  pack.  The  cell  with  the 
minimum  remaining  cell  capacity  does  not  have  the  minimum 
capacity  that  can  be  charged  in  the  majority  of  the  battery  packs,  so 
pack  capacity  is  usually  less  than  all  of  the  cell  capacities  in  the 
pack. 

For  dissipation  cell  equalization  such  as  resistor  shunt  equal¬ 
ization,  the  cell  with  minimum  cell  capacity  is  fully  charged  and 
discharged.  Thus,  pack  capacity  is  the  minimum  cell  capacity, 


Cpack(t)  =  min(C(t))  (2) 

For  non-dissipation  cell  equalization  [16,17]  which  transfers  the 
energy  from  cell  to  cell  at  all  time,  the  maximum  pack  capacity  the 
battery  pack  may  reach  is  the  mean  cell  capacity, 


Cpack(t)  =  mean(C(t))  (3) 

Therefore,  the  relationship  between  pack  capacity  and  cell 
capacities  is 


Cpack(t) 


min(Cr(t)  +  min((l  -  SOC (t)))-C(t))  without  equlization 
min(C(t))  dissipation  equlization 

mean(C(t))  non  -  dissipation  equlization 


(4) 
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No  matter  cell  equalization  is  applied  or  not,  statistical  cell 
capacities  are  needed  to  calculate  pack  capacity.  However,  cell 
capacities  are  impossible  to  be  measured  in  vehicle  applications. 
Even  though  hundreds  of  cell  capacities  were  measured  before  cell 
grouping,  we  could  just  know  initial  cell  capacities  and  initial  pack 
capacities.  As  the  cell  capacity  fade  rates  are  different  from  cell  to 
cell  due  to  the  complex  thermoelectric  conditions  in  EVs,  it  is 
almost  impossible  to  accurately  predict  every  single  cell  capacities. 

2.2.  Charging  cell  voltage  curves  of  LiFeP04 

A  typical  LiFePC>4  cell  CCVC  is  shown  in  Fig.  1.  A  cell  enters 
charging  period  after  rest  in  the  fully  discharged  state.  A  standard 
CC— CV  charging  profile  is  used  to  charge  the  cell.  The  first  period  is 
the  constant  current  charging  period  (CC  period).  The  cell  is 
charged  at  1/3C  constant  current  until  the  charge  cutoff  voltage 
which  is  3.6  V  for  LiFePC>4  cell  as  recommended.  After  CC  charging, 
the  cell  experienced  a  constant  voltage  charging  period  (CV  period) 
until  current  reaches  the  cutoff  current  which  is  1  / 20C.  The  charged 
capacity  is  divided  into  two  parts  according  to  the  charging  periods, 
CC  capacity  and  CV  capacity.  Under  the  standard  charging  profile, 
CV  capacity  accounts  for  approximately  2%  of  cell  capacity  (see 
Section  4.1 ).  As  the  CV  control  technology  of  single  cells  in  the  pack 
is  relatively  complex,  many  EVs  prefer  charging  without  CV  period. 

Fig.  2  shows  CCVCs  of  a  battery  pack  composed  of  four  LiFePC>4 
cells  in  series  with  a  nominal  capacity  of  6.5  A  h.  The  pack  starts 
charging  at  a  constant  current  of  2.16  A  which  is  1/3C  from  the  fully 
discharged  state  of  Cell  2,  and  stops  when  Cell  4  reaches  the  charge 
cutoff  voltage  of  3.6  V.  Generally,  CCVCs  are  close  to  each  other  in 
Fig.  2.  Details  of  Fig.  2  in  the  blue  dotted  and  brown  solid  boxes 
show  that  only  Cell  2  is  at  its  fully  discharge  state  before  charging 
and  only  Cell  4  is  charged  to  the  charge  cutoff  voltage.  We  can 
deduce  that  Cell  2  has  the  minimum  remaining  capacity  when 
discharged  and  Cell  4  has  the  minimum  absorbable  capacity  when 
charged.  Thus,  the  pack  capacity  is  less  than  all  of  the  4  cell 
capacities.  We  can  calculate  pack  capacity  according  to  Formula  (1 ). 
As  the  pack  starts  charging  from  the  fully  discharged  state  and  ends 
at  the  fully  charged  state,  the  pack  capacity  can  be  calculated 
directly  (For  interpretation  of  the  references  to  colour  in  this  figure 
legend,  the  reader  is  referred  to  the  web  version  of  this  article.). 

Compared  to  the  complete  charging  in  Fig.  2  which  normally 
happens  in  the  laboratory,  CCVCs  in  EVs  are  always  incomplete. 
This  is  because  drivers  always  charge  the  packs  before  they  are 
empty,  and  they  prefer  battery  packs  in  the  fully  charged  state.  It 
means  CCVCs  in  EVs  never  start  at  the  fully  discharged  state  but 
usually  end  at  the  fully  charged  state.  Fig.  3  shows  CCVCs  of  the 
same  pack  in  Fig.  2  that  simulates  possible  CCVCs  of  battery  packs 


installed  in  EVs.  The  pack  starts  charging  when  none  of  the  cells  are 
fully  discharged  and  ends  when  Cell  4  reaches  the  charge  cutoff 
voltage.  In  this  case,  the  pack  capacity  is  no  longer  measurable 
using  ampere-hour  integral. 

As  Cell  4  reaches  the  charge  cutoff  voltage  before  other  cells,  we 
manually  discharge  Cell  4  slightly.  Fig.  4  displays  CCVCs  after  cell 
equalization.  Details  of  Fig.  4  in  the  blue  dotted  and  brown  solid 
boxes  show  that  Cell  4  reaches  the  discharge  cutoff  voltage  together 
with  Cell  2  when  fully  discharged,  and  the  pack  still  stops  charging 
when  Cell  4  reaches  charge  cutoff  voltage.  This  indicates  a  cell 
capacity  difference  between  Cell  2  and  4.  We  also  infer  that  Cell  4 
has  a  minimum  cell  capacity.  Since  Cell  4  can  be  fully  charged  and 
discharged  in  the  pack,  the  pack  capacity  is  equal  to  Cell  4.  For  Cell  1 
and  3,  although  they  have  remaining  capacities  when  Cell  4  is  fully 
discharged,  they  do  not  reach  the  charge  cutoff  voltage  when  Cell  4 


Time/s  x-|o4 


Fig.  1.  A  typical  LiFeP04  cell  CCVC. 


Fig.  3.  Incomplete  CCVCs  of  the  battery  pack. 
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Fig.  4.  Complete  CCVCs  after  cell  equalization. 


is  fully  charged.  They  will  not  constrain  the  battery  pack  capacity, 
and  no  equalization  is  required  for  Cell  1  and  3  (For  interpretation 
of  the  references  to  colour  in  this  figure  legend,  the  reader  is 
referred  to  the  web  version  of  this  article.). 


3.  Uniform  CCVC  hypothesis  and  preliminary  validation 

3.1.  Uniform  CCVC  hypothesis 

The  relationship  between  OCV  and  SOC  of  LiFePC>4  is  unvaried 
[18,19]  as  OCV  is  the  external  performance  of  the  cell  electro¬ 
chemical  equilibrium  state.  For  cells  of  the  same  type,  SOC— OCV 
curves  are  uniform.  CCVC  and  OCV  curves  of  a  LiFeP04  cell  are 
compared  in  Fig.  5.  CCVC  is  very  close  to  OCV  except  at  the  end  of 


SOC 

Fig.  5.  Comparison  of  CCVC  and  OCV  curves  of  a  LiFeP04  cell. 


the  charging  process  where  CCVC  raises  to  a  high  level  due  to 
charging  polarization.  The  difference  between  charging  cell  voltage 
and  OCV  is  the  overpotential,  and  it  can  be  divided  into  three 
categories,  namely  activation  overpotential  tj3,  concentration 
overpotential  pc  and  resistance  overpotential  pr. 

According  to  Butler— Volmer  equation 

/  =  A-i0-jexpp^77aJ -expj-^^jj  (5) 

where 

/:  electrode  current 
A:  electrode  active  surface  area 
io:  exchange  current  density 
R:  universal  gas  constant 
T:  absolute  temperature 

n:  number  of  electrons  involved  in  the  electrode  reaction 

F:  Faraday  constant 

ac:  cathodic  charge  transfer  coefficient 

aa:  anodic  charge  transfer  coefficient 

p3:  activation  overpotential 

If  the  charging  environment  is  the  same  which  means  current 
and  temperature  of  the  cells  are  the  same,  and  if  SOCs  of  the  cells 
are  the  same  which  means  the  exchange  current  density  io  and 
electrodynamics  potential  £eq  (OCV)  of  the  cells  are  the  same,  the 
activation  overpotential  p3  depends  on  the  electrode  active  surface 
area  A,  cathodic  charge  transfer  coefficient  ac  and  anodic  charge 
transfer  coefficient  a3.  As  the  electrode  active  surface  area  of  cells  in 
the  same  batch  are  very  close  to  each  other  and  a3  =  ac  =  0.5  for 
most  reactions,  activation  overpotential  with  the  same  is  approxi¬ 
mately  the  same  from  cell  to  cell. 

Regarding  concentration  overpotential  pc,  according  to  equation 

%  =  §ln(CB/CE)  (6) 

where  Cb  is  bulk  concentration  and  Ce  is  electrode  concentration. 
When  cells  are  charged  with  the  same  SOC,  differences  of  the 
concentration  are  small.  Therefore  the  concentration  overpotential 
Pc  is  also  similar  from  cell  to  cell.  And  for  resistance  overpotential  rjr, 
if  the  cell  ohmic  resistances  are  the  same,  resistance  overpotentials 
pr  are  the  same. 

We  therefore  hypothesize  that  in  a  battery  pack  with  cells  of  one 
batch,  if  internal  resistances,  initial  remaining  cell  capacities  and 
total  cell  capacities  are  the  same,  CCVCs  of  the  cells  are  the  same, 
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Fig.  7.  Current  profile  and  CCVC  of  the  standard  capacity  test  for  Cell  A. 


Fig.  8.  CCVCs  of  two  LiFeP04  cells  with  the  horizontal  axis  of  ampere-hours. 


which  means  that  voltage— capacity  curves  are  completely  over- 
lapped.  We  also  deduce  that  by  CCVC  transformation,  CCVCs  of  cells 
with  different  internal  resistances,  initial  remaining  cell  capacities 
and  total  cell  capacities  can  be  overlapped  and  cell  capacities  can  be 
calculated.  Different  internal  resistances  are  eliminated  by  CCVC 
vertical  translations,  different  initial  remaining  cell  capacities  are 
compensated  by  CCVC  horizontal  translations  and  by  horizontal 
scaling,  influence  of  different  cell  capacities  to  CCVCs  are 
eliminated. 

Detail  procedures  are  demonstrated  in  Fig.  6.  We  assume  CCVC 
of  Cell  0  is  known  as  /o,  and  its  internal  resistance  ro,  initial  SOCo 
and  cell  capacity  Co  are  also  known.  CCVC  of  Cell  i  is  known  as/,  but 
its  internal  resistance  n,  initial  SOQ  and  cell  capacity  Q  are 
unknown.  Original  CCVCs  are  shown  in  Fig.  6(a).  By  CCVC  vertical 
translation  AU,  we  get  Fig.  6(b),  and  Fig.  6(c)  is  achieved  by  hori¬ 
zontal  scaling  k.  Finally,  by  horizontal  translation,  /  is  overlapped 
with/o  in  Fig.  6(d).  The  internal  resistance  n,  initial  SOC,  and  cell 
capacity  Ci  can  be  calculated  using  Formula  (7) 


(7) 


AD  >  0,  if/i  moves  downward  and  A U  <  0  if/i  moves  upward; 
k  >  1  if/i  zooms  out  and  k  <  1  if/i  zooms  in;  AAh  <  0,  if/i  moves  left 
and  AAh  >  0  if  /  moves  right. 


3.2.  Preliminary  experimental  validation 

Two  commercial  LiFeP04  cells  of  one  type  but  with  different 
capacities  of  60  A  h  (Cell  A)  and  70  A  h  (Cell  B)  are  selected  to 
preliminary  verify  the  uniform  CCVC  hypothesis.  Standard  capacity 
tests  are  used  to  test  cell  capacities,  and  also  CCVCs  are  achieved 
during  the  test.  Current  profile  and  CCVC  of  the  standard  capacity 
test  for  Cell  A  are  depicted  in  Fig.  7.  The  cell  is  charged  at  the 
constant  current  of  1/3C  (20  A)  after  it  reaches  the  discharge  cutoff 


Table  1 

Capacities  of  two  commercial  LiFeP04  cells. 

Cell  type  Cell  A  Cell  B 

Discharge  capacity  (A  h)  61.57  68.61 

CC  charge  capacity  (A  h)  60.55  67.17 

CV  charge  capacity  (A  h)  1.17  1.62 

Total  charge  capacity  (A  h)  61.72  68.79 


voltage  2.5  V.  The  charging  process  changes  from  CC  to  CV  when 
cell  voltage  reaches  charge  cutoff  voltage  3.6  V,  and  stops  when  the 
charge  cutoff  current  of  1  / 20C  (3  A)  is  reached.  After  that,  the  cell  is 
ready  for  the  standard  capacity  test.  And  the  cell  is  then  cycled  for  3 
times  at  standard  capacity  test  current  profile.  Ampere-hour  inte¬ 
gral  method  is  used  to  calculate  the  average  capacity  of  the  3  cycles 
for  the  cell  capacity. 

Cell  B  is  tested  with  the  same  method.  And  the  capacities  of  the 
two  cells  are  listed  in  Table  1.  Total  charge  capacity  is  the  sum  of  the 
CC  charge  capacity  and  CV  charge  capacity.  The  cell  capacities  are 
61.57  A  h  and  68.61  A  h  respectively  regarding  the  discharge 
capacities.  As  the  total  charge  capacities  are  close  to  the  discharge 
capacities,  the  error  is  less  than  0.2%  when  total  charge  capacities 
are  used  to  represent  the  cell  capacities. 

Fig.  8  shows  CCVCs  in  one  complete  standard  charging  cycle  of 
the  two  cells  with  the  horizontal  axis  of  ampere-hours.  The  curves 
are  not  overlapped  due  to  the  different  cell  capacities. 

We  assume  that  capacity  of  Cell  A  is  known,  and  thus  CCVC  of 
Cell  A  (/a)  is  the  standard  CCVC.  By  curve  transformation,  we  try  to 
overlap  CCVC  of  Cell  B  (fb )  with/a./a  and  the  transformed  curve  of/b 
are  shown  in  Fig.  9,  and  the  transformation  parameters  achieved  by 
trial  and  error  method  are  AU=0,k=  1.095,  AAh  —  -0.5  A  h.  Better 
parameters  are  possible,  but  the  achieved  parameters  are  good  to 
overlap  the  CCVCs.  We  use  Formula  (7)  to  calculate  the  capacity  and 
initial  SOC  of  Cell  B.  The  results  are  compared  in  Table  2.  The 
capacity  estimation  errors  are  less  than  2%  and  initial  SOC  error  less 
than  1%.  The  result  shows  that  it  is  feasible  to  estimate  cell  capac¬ 
ities  by  curve  transformation,  and  also  the  uniform  CCVC  hypoth¬ 
esis  is  reliable. 


Fig.  9.  CCVC  of  Cell  A  and  transformed  CCVC  of  Cell 
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Table  2 

Cell  capacity  and  initial  SOC  estimation  errors  by  curve  transformation. 

Cell  type  Cell  A  Cell  B  Cell  B  Difference/ 

(estimated)  (measured)  accuracy 
CC  charge  60.55  66.30  67.17  0.87/1.30% 

capacity  (A  h) 

Total  charge  61.72  67.58  68.79  1.21/1.76% 

capacity  (A  h) 

Initial  SOC  0  -0.8%  0  0.8% 


The  main  purpose  of  this  paper  is  to  estimate  the  pack  capacity 
and  cell  capacity,  so  internal  resistances  are  not  concerned  in  this 
study.  It  is  better  to  reduce  the  three-dimensional  problem  (All, 
AAh  and  k)  to  two-dimensional  problem  (AAh  and  k).  The  slope  of 
CCVC  is  not  influenced  by  CCVC  vertical  translation,  and  it  is 
independent  of  parameter  All.  We  use  the  parameters  AAh  and  k  to 
overlap  the  slope  curve  of  CCVC  which  is  called  as  voltage— capacity 
rate  curve  (VCRC)  in  this  paper  for  convenience.  The  voltage- 
capacity  rate  stands  for  the  voltage  increase  in  an  ampere-hour. 
Fig.  10  shows  VCRCs  of  cells,  the  CCVCs  of  which  are  plotted  in 
Fig.  8.  As  expected,  VCRCs  are  similar.  With  horizontal  translation 
AAh,  horizontal  contraction  k  and  vertical  expansion  k,  VCRCs  will 
overlap. 

As  we  fully  charge  and  discharge  the  cells,  the  initial  SOCs  of  the 
cells  are  the  same.  So  we  do  not  need  to  focus  on  the  initial  SOC  in 
our  preliminary  experimental  validation  and  no  horizontal  trans¬ 
lation  is  required.  Horizontal  contraction  k  and  vertical  expansion  k 
are  applied  to  transform  VCRC  of  Cell  B.  Fig.  11  demonstrates 
comparison  of  the  transformed  VCRCs  of  Cell  B  at  different 
contraction/expansion  k  with  VCRC  of  Cell  A.  when  k  =  1.1,  the 
transformed  VCRC  of  Cell  B  is  very  close  to  VCRC  of  Cell  A.  The  blue 
dashed  box  shows  that  when  k  —  1.095  or  k  =  1.105,  the  trans¬ 
formed  VCRCs  of  Cell  B  are  not  as  good  as  when  k  =  1.1.  We  conclude 
that  k  -  1.1  is  a  better  solution,  and  the  corresponding  calculated  CC 
charge  capacity  of  Cell  B  is  66.61  A  h,  which  is  0.56  A  h  less  than  the 
measured  capacity.  The  relative  error  is  less  than  1%  (For  inter¬ 
pretation  of  the  references  to  colour  in  this  figure  legend,  the  reader 
is  referred  to  the  web  version  of  this  article.). 

The  accuracy  of  the  cell  capacity  estimation  by  CCVC/VCRC 
transformation  is  proved  to  be  satisfactory.  However,  it  is  subjective 
by  observing  the  overlapped  CCVC/VCRC  and  the  trial  and  error 
method  also  causes  a  lot  of  time,  so  it  cannot  be  applied  in  EVs  for 
cell  capacity  estimation.  In  addition,  CCVCs  of  the  two  cells  are 
achieved  respectively  and  completely,  but  CCVCs  in  EVs  are  usually 
not  complete  and  errors  of  the  estimated  cell  capacities  and  pack 
capacity  are  still  to  be  verified. 


Ah /Ah 


Fig.  11.  VCRC  of  Cell  A  and  the  transformed  VCRCs  of  Cell  B  with  different  trans¬ 
formation  parameter  k. 


4.  Pack  capacity  estimation 

4.1.  Charging  experiment  of  a  battery  pack 

A  small  battery  pack  consists  of  four  LiFePCU  cells  in  series  with 
a  nominal  capacity  of  6.5  A  h  is  employed  to  simulate  possible 
CCVCs  of  battery  packs  installed  in  EVs.  Cell  capacities  are  tested 
before  grouped,  and  one  of  the  cell  capacity  test  profile  is  shown  in 
Fig.  1.  The  cell  capacities  are  listed  in  Table  3.  The  maximum 
capacity  of  the  four  cells  is  Cell  3,  and  the  minimum  capacity  Cell  4. 
The  CC  charge  capacity  is  approximately  98%  of  the  total  charge 
capacity.  The  cells  are  at  the  fully  charged  state  after  the  standard 
capacity  test.  But  in  order  to  simulate  different  SOCs  in  the  pack,  the 
initial  remaining  capacities  are  conditioned  at  a  discharge  current 
of  3  A.  The  cells  are  discharged  for  5,  6,  9  and  0  min  with  0.25,  0.3, 
0.45  and  0  A  h  released  from  the  respective  cells.  The  condition  can 
be  arbitrary  as  SOCs  in  the  pack  are  various.  The  conditioned  initial 
remaining  capacities  are  also  listed  in  Table  3. 

After  cell  capacity  test  and  conditioning,  the  cells  are  connected 
in  series  for  a  complete  pack  charging  in  Fig.  2.  As  the  initial 
remaining  capacity  of  Cell  2  is  the  minimum,  it  reaches  the 
discharge  cutoff  voltage  before  other  cells.  Cell  4  reaches  the  charge 
cutoff  voltage  before  others  as  it  has  the  minimum  absorbable 
capacity.  We  can  directly  calculate  the  pack  capacity  by  ampere- 
hour  method  as  the  pack  has  a  complete  charge.  The  pack 
capacity  is  7.033  A  h,  which  is  less  than  any  cell  capacities  as 
expected. 

As  CCVCs  in  EVs  never  start  at  the  fully  discharged  state  but 
usually  end  at  the  fully  charged  state,  we  further  simulate  possible 
CCVCs  of  battery  packs  installed  in  EVs  in  Fig.  3.  The  pack  starts 
charging  when  none  of  the  cells  are  fully  discharged  and  ends  when 
Cell  4  reaches  the  charge  cutoff  voltage.  The  charged  ampere-hours 
of  6.226  A  h  is  no  longer  the  pack  capacity. 

Due  to  the  unknown  cell  capacities,  cell  SOC  estimation  by 
ampere-hour  method  is  inaccurate,  but  cell  SOC  estimation  by 
SOC— OCV  curve  does  not  depend  on  cell  capacities.  The  SOC— OCV 
curve  of  the  LiFeP04  is  plotted  in  Fig.  12.  Probability  density 
distribution  of  OCVs  is  drawn  along  the  vertical  axis.  As  can  be  seen, 


Table  3 

Cell  capacities  and  initial  remaining  capacity  of  4  LiFePCU  cells. 

Cell  Cell  #1  Cell  #2  Cell  #3  Cell  #4 

Total  charge  capacity  (A  h)  7.6085  7.4494  7.6684  7.2831 

CC  charge  capacity  (A  h)  7.4464  7.2934  7.5191  7.1171 

Initial  remaining  capacity  (A  h)  7.3585  7.1494  7.2184  7.2831 


Fig.  10.  Slope  i 


of  CCVCs  of  the  2  LiFeP04  cells  (VCRCs). 
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Fig.  12.  SOC— OCV  curve  of  the  LiFeP04  Cell  and  the  probability  density  distribution  of 
OCVs. 

when  SOC  is  less  than  15%  and  the  voltage  is  less  3.214  V,  the  OCV 
probability  density  is  small  which  means  that  the  accuracy  of  the 
estimated  SOC  is  high  when  the  measured  OCV  is  less  than  3.214  V. 
But  if  the  measured  OCV  is  larger  than  3.214  V,  the  OCV  probability 
density  is  large  and  its  corresponding  SOC  interval  is  large  which 
means  the  accuracy  of  the  estimated  SOC  is  not  reliable. 


Ah 


Fig.  13.  Objective  function  construction  for  VCRC  transformation. 


The  cell  capacity  of  VCRCgo(Ah)  is  already  known.  The  objective 
is  to  determine  the  best  transformation  parameter  T  =  [k,  AAh]  to 
make  VCRC  transformation  from  VCRC  fi(Ah)  to  VCRC 
g'j(Ah)  =  kgi(Ah/k  -  AAh)  so  as  to  overlap  VCRC  g'(Ah)  and  VCRC 
go(Ah).  As  demonstrated  in  Fig.  13,  an  approximate  method  for  the 
best  curve  overlapping  is  to  minimize  the  distances  of  the  corre¬ 
sponding  points  on  the  curves.  The  coincident  domain  of  go(Ah) 
and  g[(Ah)  is  denoted  as  D  =  [Ahmin  Ahmax],  Corresponding  points 
are  chosen  every  certain  ampere-hour,  e.g.  (Ahmax  -  Ahmin)/50. 
These  points  are  marked  as  follows: 


Po  =  {Pot,  Po2,  -Pon}  =  {(Ah,,  g0(Ah, )),  (Ah2.  g0(Ah2)),  ...(Ah„,  g0(Ah„))}; 
P'i  =  {p[v  Pw-P'in}  =  UAhi>  g;(Ah1)),(Ah2,  gf(Ah2)),...(Ah„,  gj(Ahn)) } ; 
Ah,,  Ah2...Ah„,eD 


Before  the  incomplete  charging  shown  in  Fig.  3,  the  battery  pack 
rests  for  3  h  to  achieve  cell  OCVs.  All  the  measured  OCVs  are  less 
than  3.21  V,  and  accurate  SOCs  can  be  achieved  by  SOC-OCV  curve 
in  Table  4.  As  Cell  4  reaches  the  charge  cutoff  voltage,  the  SOC  is 
100%  when  the  pack  is  fully  charged.  Capacity  of  Cell  4  then  can  be 
calculated:  6.226  A  h/(100%— 13.2%)  =  7.17  A  h.  Compared  to  the 
measured  capacity,  the  calculated  capacity  of  Cell  4  is  0.05  A  h 
larger.  The  error  is  less  than  1%. 

However,  only  capacity  of  Cell  4  can  be  calculated,  because  it  is 
the  only  cell  that  is  fully  charged.  For  other  three  cells,  SOCs  are 
indeterminate  at  the  end  of  charging,  so  the  capacities  cannot  be 
calculated  as  Cell  4. 

The  estimated  capacity  of  Cell  4  is  the  foundation  for  capacity 
estimation  of  the  other  cells.  With  the  known  capacity  of  Cell  4, 
CCVC  transformation  can  be  applied  to  the  other  cells.  Capacities  of 
these  cells  as  well  as  the  pack  capacity  can  be  estimated. 

4.2.  Pack  capacity  estimation  by  genetic  algorithm 

We  have  verified  that  VCRC  transformation  is  reliable  for  cell 
capacity  estimation.  But  we  need  to  seek  a  better  mathematical 
solution  as  the  trial  and  error  method  is  ineffectiveness. 


The  objective  function  is  constructed  in  Eq.  (8) 

G=  £|gi(Ahj)— g0(Ah,)|  (8) 

Theoretically  when  the  best  transformation  parameter  T  — 
[k,  AAh]  is  applied  to  VCRC  g,,  g[  is  overlapped  with  VCRC  go.  As 
a  result,  Po  and  Pj  are  the  same,  so  the  objective  function 
G  =  Gmin,giobai  —  0  is  the  global  minimum.  As  the  voltage 
measurement  error  always  exists,  there  is  no  doubt  that  g[  cannot 
completely  match  VCRC  go  no  matter  what  transformation 
parameter  T  is  applied.  But  if  the  objective  function  G  reaches  its 
minimum  value,  g[  matches  VCRC  g0  in  the  maximal  degree.  The 
transformation  parameter  T  is  then  considered  suitable  for  VCRC 
transformation.  We  can  get  the  cell  capacity  and  its  SOC  corre¬ 
sponding  to  VCRC  g,  using  Formula  (7)  with  T  =  [k,  AAh] 
substituted. 

The  transformation  parameter  T  =  [k,  AAh]  can  be  determined 
by  some  global  optimization  methods.  Because  the  objective 
function  G  is  not  continuously  differentiable,  it  is  difficult  to  utilize 
the  traditional  optimization  methods.  Though  the  exhaustive 
method  is  a  possible  approach,  it  is  lack  of  efficiency.  Holland  and 


Table  4 

Cell  SOCs  before  charging  and  the  determined  cell  capacity. 

Cell  Cell  #1  Cell  #2  Cell  #3  Cell  #4  Pack 

Measured  CC  charge  7.4464  7.2934  7.5191  7.1171  7.033 

capacity  (A  h) 

Initial  SOC  (it)  13.6  11.1  11.7  13.2 

Calculated  CC  charge  —  —  —  7.17  — 

capacity  (A  h) 


Table  5 

Settings  in  Matlab®  GA  toolbox. 
Parameter 

Number  of  variables 

Population  size 
Number  of  generations 
Selection  function 
Crossover  fraction 
Mutation  function 


Value 

Formula  (8) 

2 

ke  [0.8  1.2],  AAhe  [— 1  1] 


Stochastic  uniform 
0.8 

Adaptive  feasible 
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Table  6  Table  7 

Estimated  cell  capacities  and  SOCs  by  GA  and  their  errors.  The  estimated  pack  capacity  and  its  error. 


Cell  #2 


Cell  #3 


K 

Calculated  CC  charge  capacity  (A  h) 
Measured  CC  charge  capacity  (A  h) 
Capacity  estimation  error  (%) 
AAh(Ah) 

Calculated  SOC  (%) 

Measured  SOC  (%) 

SOC  estimation  error  (%) 


1.044  1.021  1.052 

7.49  732  7.54 

7.4464  7.2934  7.5191 

0.59  0.36  0.28 

-0.011  -0.124  -0.090 

13.06  11.47  11.95 

13.6  11.1  11.7 

0.54  0.37  0.25 


Cr(Ah) 

(1  —  SOC)*C  (A  h) 
min(C,( t))  (A  h) 
min((l  -  SOC(t)).C(t))  (A  h) 
Calculated  Cpaci(  (A  h) 
Measured  CPack  (A  h) 

Cpack  estimation  error  (%) 


0.98  0.84  0.90 

6.51  6.48  6.64 

0.84 

6.22 

7.06 

7.033 

0.4 


Bold  value  indicates  they  are  the  smallest  numbers  of  each  rows. 


0.95 

6.22 


his  students  in  the  University  of  Michigan,  USA,  proposed  the 
Genetic  Algorithm  (GA)  which  is  an  efficient  and  strong  robust 
global  optimization  algorithm.  It  can  quickly  and  reliably  solve 
problems  with  complex  objective  functions  which  are  not  contin¬ 
uously  differentiable  [20],  GA  is  employed  to  search  the  global 
optimization  transformation  parameter  T  in  this  paper. 

GA  is  based  on  natural  selection  which  is  completely  different 
from  the  traditional  optimization  methods.  GA  simulates  natural 
biological  evolution  according  to  the  “fitness  level”  of  the  individ¬ 
uals  which  are  a  large  set  of  possible  solutions  to  a  given  problem. 
Through  the  genetic  operators  (selection,  crossover,  and  mutation) 
and  natural  selection,  better  generations  are  bred.  By  selection 
operator,  parent  solutions  which  have  better  fitness  level  are  more 
likely  to  reproduce  which  means  better  genes  are  more  likely  to 
dominate  the  next  generation.  Crossover  combines  the  features  of 
two  parents  at  a  certain  crossover  fraction  to  form  new  solutions  by 
swapping  corresponding  segments  of  parent  chromosomes.  By 
randomly  changing  one  or  more  genes  at  a  low  mutation  rate, 
mutation  introduces  variability  into  the  next  generation  that  will 
stop  GA  converging  at  a  local  minimum. 

Matlab®  GA  toolbox  is  employed  to  resolve  the  global  optimi¬ 
zation  VCRC  transformation  parameter  T  in  this  paper.  The 
parameters  for  the  VCRC  transformation  are  listed  in  Table  5.  The 
fitness  function  is  the  objective  function  in  Eq.  (8).  In  the  Matlab® 
GA  toolbox,  a  smaller  fitness  function  value  means  a  better  fitness 
level.  As  a  result,  the  best  individual  of  all  the  generations  is  the 
best  transformation  parameter  to  minimize  the  objective  function 
G.  We  have  2  variables  in  the  transformation  parameter  T.  Their 
bounds  are  set  fee  [0.8  1.2],  AAhe [-1 1],  because  the  differences  of 
cell  capacities  and  SOCs  will  not  exceed  20%  normally  when  there 
are  no  faults  occurred  in  the  pack.  100  generations  are  bred  and 


the  population  size  is  100  for  each  generation.  Selection  function 
uses  stochastic  uniform  algorithm  which  ensures  that  the  fitness 
levels  of  the  parent  individuals  are  corresponding  to  their  selec¬ 
tion  probability.  Mutation  function  uses  adaptive  feasible  algo¬ 
rithm  which  generates  mutation  direction  randomly  and  selects 
the  mutation  step  length  in  the  range  of  the  transformation 
parameter  bounds. 

Transformation  parameters  of  VCRCs  for  Cell  1,  2  and  3  are 
optimized  by  GA  and  the  results  are  shown  in  Table  6.  The  capacity 
and  SOC  estimation  errors  are  all  well  below  1%  which  indicates 
that  GA  is  considerably  effective  for  parameter  optimization  of 
VCRC  transformation. 

Take  Cell  1  for  example,  the  transformation  parameter 
T  =  [1.044  -0.011],  i.e.  k  =  1.044,  AAh  =  -0.011.  According  to 
Formula  (7),  the  estimated  cell  capacity  is  7.49  A  h  with  an  error  of 
0.59%  and  the  estimated  SOC  is  13.06%  with  an  error  of  0.54%. 
Fig.  14  depicts  VCRC  and  CCVC  of  Cell  1  before  and  after  curve 
transformation.  VCRC  and  CCVC  of  Cell  4  are  also  plotted  in  Fig.  14 
using  the  red  solid  curves.  The  blue  dashed  curves  are  VCRC  and 
CCVC  before  transformation,  and  after  transformation  they  turn  to 
blue  solid  curves.  As  can  be  seen,  the  transformed  VCRC  and  CCVC 
of  Cell  1  are  well  overlapped  with  those  of  Cell  4  in  the  figure  (For 
interpretation  of  the  references  to  colour  in  this  figure  legend,  the 
reader  is  referred  to  the  web  version  of  this  article.). 

Formula  (1)  is  applied  to  calculate  the  pack  capacity  based  on 
the  results  of  cell  capacities  and  SOCs.  The  estimated  pack  capacity 
is  shown  in  Table  7.  As  expected  in  Section  2.2,  Cell  2  has  the  least 
ampere-hours  that  can  be  released  and  Cell  4  has  the  least  ampere- 
hours  that  can  be  charged.  The  estimated  pack  capacity  is  7.06  A  h 
with  an  error  of  0.4%  compared  to  the  measured  pack  capacity 
7.033  Ah. 


a  b 


Fig.  14.  The  comparison  of  VCRC  and  CCVC  of  Cell  4  and  those  of  Cell  1  before  and  after  i 


transformation. 
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5.  Conclusion 

An  accurate  estimation  of  the  battery  pack  capacity  is  significant 
for  forecasting  the  EV  driving  range.  But  because  of  the  diversiform 
driving  conditions  and  the  cell  variations,  it  is  difficult  to  accurately 
determine  battery  pack  capacities  in  EVs  by  model  prediction  or 
direct  measurement.  This  paper  studies  CCVC  and  VCRC  trans¬ 
formation  for  the  estimation  of  the  LiFeP04  battery  pack  capacities 
in  EVs. 

Firstly,  we  believe  that  with  cell  variations  taking  into  consid¬ 
eration  the  pack  capacity  is  different  from  the  cell  capacities  in  the 
pack.  Thus  we  deduce  the  quantitative  relationship  between  pack 
capacities  and  cell  capacities  in  the  packs. 

Secondly,  we  propose  the  uniform  CCVC  hypothesis  that  in 
a  battery  pack  with  cells  of  one  batch,  if  internal  resistances,  initial 
remaining  cell  capacities  and  total  cell  capacities  are  the  same, 
CCVCs  of  the  cells  are  the  same.  We  estimate  cell  capacities  by 
overlapping  the  CCVCs  using  CCVC  transformation.  CCVCs  of  two 
LiFePC>4  cells  with  large  capacity  difference  are  used  to  verify  the 
hypothesis  and  the  result  shows  a  satisfactory  accuracy  of  cell 
estimation. 

Finally,  we  further  develop  an  equivalent  simplified  approach 
using  VCRC  and  implement  GA  to  find  the  optimum  transformation 
parameter  for  overlapping  VCRCs.  A  small  battery  pack  with  four 
LiFeP04  cells  in  series  is  employed  to  verify  the  method  and  the 
result  shows  that  the  estimation  errors  of  both  pack  capacity  and 
cell  capacities  are  less  than  1%. 

With  the  proposed  method  and  easily  achieved  CCVCs  in  EVs, 
the  pack  capacity  can  be  precisely  estimated.  In  addition,  cell 
capacities  calculated  in  battery  packs  will  support  the  study  of  cell 
degradation  and  variations  in  vehicle  driving  conditions,  and  the 
calculated  cell  capacities  are  also  the  key  information  for  cell 
equalization  strategies. 

Further  studies  are  to  be  carried  on  in  the  following  three 
aspects.  First,  more  CCVCs  in  EVs  are  encouraged  to  prove  the 
uniform  CCVC  hypothesis  as  the  authors  have  limited  data;  Second, 
cells  of  other  types,  e.g.  LiMn04,  are  to  be  studied  besides  LiFePC>4 
cells  to  discover  whether  the  uniform  CCVC  hypothesis  is  suitable 
for  other  cells;  Third,  more  simple  algorithms  are  to  be  developed 


for  on  line  estimation  of  cell  capacities  and  pack  capacities  in  EVs  as 
GA  is  still  complex  in  computing  for  vehicle  applications. 
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